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The two-flavor color superconductivity is studied over a wide range of baryon density with a 
O . single model. We pay a special attention to the spatial-momentum dependence of the gap and to 

the spatial-structure of Cooper pairs. At extremely high baryon density (~ O(10 10 po) with po being 



in 



the normal nuclear matter density), our model becomes equivalent to the usual perturbative QCD 
treatment and the gap is shown to have a sharp peak near the Fermi surface due to the weak-coupling 
| nature of QCD. On the other hand, the gap is a smooth function of the momentum at lower densities 

(~ O(10po)) due to strong color magnetic and electric interactions. To study the structural change 
£N| ' of Cooper pairs from high density to lower density, quark correlation in the color superconductor 

^ \ is studied both in the momentum space and in the coordinate space. The size of the Cooper pair 

is shown to become comparable to the averaged inter-quark distance at low densities. Also, effects 
of the momentum-dependent running coupling and the antiquark pairing, which are both small at 
high density, are shown to be non-negligible at low densities. These features are highly contrasted 
to the standard BCS superconductivity in metals. 
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I. INTRODUCTION 



45 
D., 

Because of the asymptotic freedom and the Debye screening in QCD, deconfmed quark matter is expected to be 
realized for baryon densities much larger than the normal nuclear matter density Q. Furthermore, any attractive 
quark-quark interaction in the cold quark matter causes an instability of the Fermi surface by the formation of Cooper 
pairs and leads to the color superconducting phase [|| |6| . 

Current understanding on the color superconductivity has been based on two different theoretical approaches. One 
is an analysis by the one-loop Schwinger-Dyson equation with perturbative one gluon exchange, which is valid at 
asymptotically high densities ^] ■ The dominant contribution to the superconducting gap at very high density 
comes from the collinear scattering through the long range magnetic gluon exchange [p_0| . In such a weak-coupling 
regime, formation of Cooper pairs takes place only in a small region near the Fermi surface. The other approach to 
the color superconductivity is the mean-held approximation with QCD inspired 4-Fermi model which is introduced 
to study lower density regions |^|,|5|,[ll],[l2| . In this approach, magnitude of the gap becomes as large as 100 MeV and 
is almost constant in the vicinity of Fermi surface as a function of momentum. 

The main purpose of this paper is to calculate the superconducting gap over a wide range of baryon density with 
a single model and to make a bridge between high and low density regimes. In particular, we make an extensive 
analysis on the structural change in spatial- momentum dependence of Cooper pair from high to low densities. At 
relatively low baryon densities, the gluonic attraction becomes effective for all quarks inside the Fermi sea and sizable 
Cooper pairing takes place for wide range of momentum away from the Fermi surface. This is in contrast to the 
phonon interaction with an intrinsic Debye cutoff in the BCS-type superconductivity Jl3[ . The spatial-vaovaeatxym 
dependence of the gap, diffuseness of the Fermi surface, quark-quark correlations in the superconductor and the spatial 
size of Cooper pairs are the characteristic quantities which reflect the departure from the weak-coupling picture. So 
far, only a few investigations have been made along this line P, |l4| , |l5f . 

Throughout the paper, we limit our discussions to 2-fiavor color superconductivity partly because our primary 
interests are in low densities and partly because the analysis is simpler than in the three-flavor case fl6[| . To investigate 
how the weak-coupling picture is modified at low densities in a qualitative manner, we adopt a model called the 
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improved ladder approach in which the one-loop Schwinger-Dyson equation with infrared safe running coupling is 
used. This model is known to reproduce the physics of the QCD vacuum and the qq meson properties reasonably well 
[ L7-p~9[| . For recent applications of similar kind of models to QCD at finite temperature and density, see the review 
[ 21 1 . As we will show later, the model can naturally reproduce the correct asymptotic behavior of the superconducting 
gap, in the high density limit. Therefore, it is a suitable model for our purpose to study the color superconductivity 
for a wide range of density. To make our analysis as simple as possible, we use a gluon propagator with a minimal 
static and dynamical corrections (Debye screening and Landau damping). The similar form of the gluon propagator 
is also adopted in Ref. 22 1. 

This paper is organized as follows. In Sect. II, we define our model and derive relevant gap equations. We first 
consider the weak-coupling region at high density and discuss the general properties of the gap as a function of the 
spatial- momentum. In Sect. Ill, we solve the momentum-dependent gap equation numerically for a wide range of 
density. A large structural change of the momentum-dependent gap will be shown. The quark and antiquark occupa- 
tion numbers, the correlation of quarks in the color superconductor, and the coherence length are also investigated. 
Summary and concluding remarks are given in the last section. In Appendix A, the gauge dependence of the gap 
equation is discussed. In Appendix B, we derive the asymptotic behavior of the spatial quark-quark correlation in 
the weak coupling. 



II. GAP EQUATION WITH SPATIAL MOMENTUM DEPENDENCE 

In this section, we define our model and derive relevant gap equations in two flavor QCD with the one-gluon exchange 
interaction. We start with the standard Schwinger-Dyson equation for the quark self-energy with four momentum 
dependence 0. Then it is reduced to a gap equation with spatial-momentum dependence. General properties of 
the solution of our gap equation in the weak-coupling limit are also examined. The gap with spatial-momentum 
dependence enables us to study the diffuscncss of the Fermi surface and the the structure of Cooper pairs in the color 
superconductor. 

Throughout this paper, we use the following notation for the four and three momenta : k p = (k°, k), k = \k\, k = 
k/k. Also, we limit ourselves to the system at zero temperature and work in the Minkowski space with the quark and 
gluon propagators obeying the Feynman boundary condition. 



A. The model 



Let us start with the definition of quark self-energy and the superconducting gap. Using the standard Nambu- 
Gor'kov formalism with a two component Dirac spinor W = (-0,-0*), the quark self-energy S(fc p ) with the Minkowski 
4-momentum k p is written as 

E(fe p ) = 5 - 1 (fc p )-S , - 1 (M=( A( ° M A( o M ). (2- 1 ) 

where the superconducting gap A and A = 7oA'l'7o enter through the off-diagonal components of E. The diagonal 
components of £ is neglected throughout this paper 0. S(k p ) and So(k p ) are the full and free quark propagators, 
respectively. The free propagator is taken to be a form SQ 1 (k p ) = diag.(^ + /I, — f} T ). We ignore the mass of u 
and d quarks. 

In the ladder approximation of the one gluon exchange, the Schwinger-Dyson equation for £ is written as 

S(fcp) = -ij 70 a 2 KS{q p )TlDf v {q P ~ M> (2-2) 

where is the gluon propagator in medium and is the quark-gluon vertex. 

To study the flavor anti-symmetric, color anti-symmetric and J = + channel (which is the most attractive channel 
within the one gluon exchange model), we assume the following structure of the gap function [^||J2^], 

A(k p ) = (A 2 r 2 C 75 ) (A+(fc p )A + (fe) + A_(fc„)A_(fc)) . (2.3) 

Here t 2 is the Pauli matrix acting on the flavor space, A 2 is a color anti-symmetric Gell-Mann matrix, and C is the 
charge conjugation. A±(fe) = (1 ± k ■ ot)/2 is the projector on positive (+) and negative (— ) energy quarks. A± will 
be identified as the gap in quark and antiquark channel, respectively. 



For the spinor and color structures of the vertex in Eq. (2.2), we use the bare vertex T° = diag.(7 M A a /2, —( r y fl X a /2) t ) 



§. For g 2 in Eq. we use a momentum dependent coupling g 2 (q,k) which does not diverge at low momentum 



scale. In the "improved ladder approximation" |17|, g (q, k) is taken to be 



Sir 2 



A) ln((pi ax +^)/ A 2)' 



Pn 



= max((j, k), 



(2.4) 



where /3q = (IIN C — 2iV/)/3, p 2 plays a role of a phenomenological infrared regulator, and A dictates the logarithmic 
decreases of the vertex at high momentum. The quark propagator in the improved ladder approximation in the vacuum 
is known to have a high momentum behavior consistent with that expected from the renormalization group and the 
operator product expansion. For the numerical values of A and p 2 , we adopt 400 MeV and 1.5 A 2 respectively. They 
are determined to reproduce the low energy meson properties for Nf = 2 by solving the Schwinger-Dyson equation 
for the quark propagator and the Bethe-Salpter equation for the qq bound state in the vacuum [p.8| . 
We use the following gluon propagator in medium in Eq. (2.2) in the Landau gauge, 



D pv (kp) 
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ph 



k 2 +iM 2 \k Q \/\k\ k 2 



(2.5) 



where P^ L are the transverse and longitudinal projectors satisfying the relations, P? 



kikj, Pqq 



P T 

Oi 



and 



pL 

pv 



~9ij,u + k^kv/k 2 — Pj v . The longitudinal part of the propagator has static screening by the Debye mass 
rrip = (Nf /2-K 2 )g 2 /i 2 , while the transverse part has dynamical screening by the Landau damping M 2 — (-K/^m 2 ^ [ p0| . 
The above form of the propagator is a quasi-static approximation of the full gluon propagator in the sense that only 
the leading frequency dependence is considered. The same form is also adopted in Ref. p2| . 

Here we make a brief comment on the gauge-parameter £ dependence of the gap equation. Although physical 
quantities should be gauge invariant, special truncation scheme of the diagrams such as the ladder approximation 
introduces gauge-parameter dependence. It has been claimed that (i) A+ is gauge invariant at extremely high density 
[0, and (ii) the gauge dependent contribution only begins to decrease at extraordinarily high densities \i > 10 8 MeV 
and seems to converge to the result of the Landau gauge (£ = 0) |24|]. We will show in Appendix A that, if one 
appropriately treats the momentum dependence, the ^-dependence survives in the gap equation for A + but gives a 
sub-leading contribution. Under this remark, we will consider only the Landau gauge (£ = 0) throughout this paper. 
The gap equation under the approximations shown above is obtained from the 2-1 element of the Schwinger-Dyson 



(matrix) equation (2.2). Then the final form of the gap equation for A±(k p ) — A±(ko, k) reads 



i±W~* / (04 5 2 (9,fc){^tr [k ± {k)Y^{q)Y 



A+(<Z P ) 



q 2 -E 2 + {q)-\A + {q p )\ 2 



A-(? P ) 



q 2 -E 2 _(q)-\A.(q p )\ 2 



D pv {qp k p ), 



(2.6) 



where E±(q)(= q =p /i) is the reduced energy for free quarks (for + sign) and for antiquarks (for — sign) 



B. Analytic properties of the gap in the weak coupling 



1. Gap equation in the weak-coupling limit 



Let us first study the analytic properties of the gap equation with spatial-momentum dependence. To obtain 
analytically tractable equations, we consider here the weak-coupling (high density) limit. This "weak-coupling ap- 
proximation" implies the approximation where we take only the color-magnetic interaction and the real part of A + . 
Momentum dependence of g 2 (q,k) is neglected. Contributions from the a ntiqu arks proportional to A_ and the 
momentum-dependent factor originating from the projection operator in Eq. (2^) are also neglec ted. We also ignore 
the imaginary part of the propagator because it is factor g 2 smaller than the real part. In Eq. (|]^), the frequency 
integral is performed by picking up quasi-particle poles in the quark propagator. After integrating over the angular 
variable, we obtain the following gap equation with the spatial-momentum dependence 



A+(fc) 



<r 



1447T 2 



dq 



A+(?) 



E 2 + (q) + A 2 + (q) 



(g + fc) 6 + M 4 ( e +-e+) 2 
(q - kf + M 4 (e+ - e+) 2 



+ (4 



-4) 



(2.7) 



where E + (q) = q — ji and we have defined A + (q) = A + (e+ , q) for notational simplicity, e+ is the quasi-particle energy 
as a solution of (e^~) 2 = E 2 _(q) + A^e^", q). This is the final form of an approximate gap equation in the weak-coupling 
limit. A + (fc) corresponds to the gap of the on-shell quasi-particles in the color superconductor since ko is taken to be 



2. Gap near the Fermi surface in weak coupling limit 



Let us consider the behavior of the gap near the Fermi surface. As we claimed in Introduction, the gap equation 



( p.7|) reproduces the correct high density solution near the Fermi surface, which was previously derived in Refs. [|jJlC|] 



for frequency dependent gap, and in Ref. || for momentum dependent gap. For this purpose, we pick up the most 
singular part of the integrand by taking q = k = [X in the bracket except for the term (e+ — e^) 2 in the denominator. 



Then Eq. (2.7) reduces to 



a + (* ~ /*) = 4i f * / A+(g) ln r#^< ( 2 - 8 ) 

72 ^"- /n lE\{q) + |A + (q)| 2 14 -4 I 



where b = 2(2/i/M) 2 = (2 / N f)(Z2ir / g 2 ) . This equation makes sense only in the vicinity of the Fermi surface \k— fi\ <C S 
and coincides to the one derived in Ref. 0. 

In the weak-coupling regime (g/(3\/27r) <C 1), the approximate solution of this equation is known to be |0j 

A+(fc ~ n) a A + ( M ) cos ( -2=- In |fc ~ Ml + , A+( M )=26 Aie -^ 2 /^)/f, (2.9) 



3^277 A+(/i) 

where e^T = (|fc — /i| 2 + A 2 ^^)) 1 / 2 . This solution is a symmetric function of (k — fi) and decreases as \k — /z| increases. 
Inclusion of the color-electric interaction only modifies the prefactor as b — > (2/7V/) 5 / 2 2567r 4 /g 5 , which enhances the 

gap in the weak coupling. The characteristic form of the gap e~ c l g in Eq. J2.9| ), which is not the BCS form e~ c / g2 , 
was first derived in Ref. |l0||. 



3. Gap far from the Fermi surface in weak coupling limit 



Next, we examine the behavior of the gap far from the Fermi surface. When k 3> fJ., the gap equation (2.7) can be 
cast into the following differential equation 



k ^ A +^ + 4k A ^ + &l A +^=°> 



where we have assumed that A+(fc) — > as k — > oo. The asymptotic solution of this equation for g 2 < 3ir 2 is 

A + (fc > fi) oc k~ x , A=l+ v/l -,g 2 /37r 2 . (2.11) 

Therefore, in the weak coupling g <C 1, the gap decreases slowly as 1/A; 2 for k 3> [i and does not vanish for arbitrary 
large k. This long tail is a direct consequence of the gluonic interaction which allows scattering of quarks with an 
arbitrary energy-momentum transfer. Otherwise, the gap cannot take the nonzero value far away from the Fermi 
surface. 

This is in sharp contrast to the gap in the standard weak-coupling BCS superconductivity in metals, [|] where 



1 This is a solution of the simplest gap equation discussed in the original BCS paper Jl^] 

A Bcs m= l Vy A BCS (fc / ) 

2 V ^ V(^-») 2 + ABCS (k') 2 

where Vi,y 7^ for cf — ud < £k,k' < £ p + and v = 0, otherwise. 



A BCS (p) = (n IV*- I<WD (2-12) 
v 10 otherwise, 

with e p , £f and ud being the electron energy, the electron Fermi energy and the ultraviolet Debye cutoff of the phon on 



interaction originating from the lattice structure, respectively. As we mentioned above, the long tail in Eq. (2.11) is 
intimately related to the absence of the intrinsic ultraviolet-cutoff such as cl>d in QCD. 



C. Full gap equation 



Now let us derive the full gap equation with momentum dependence from the original gap equation (2.6) with 4 
momentum dependence. 

In Eq. (|2.6| ), we adopt the static gluon propagator with the Debye screening (for the electric part) and the Landau 
damping (for the magnetic part). The imaginary part of the gap function will be neglected as we stated before. Using 
the approximation Ph, w — S^qSvq and performing the go and angular integration, the gap equation becomes 
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A±(fc) = / dq V ± (q, 
Jo 



6 q ' e k)~ 



A+(<?) 



pOO 

dqV T (q,k;e~,ef)- 



A-(«) 



VE + (q) 2 + \A + (q)\ 2 Jo + w ' 9 ' k '^ElW 

where we used a simplified notation: A±(e fc , k) —> A±(k), and V± is defined by 

g 2 (q,k)q\l ((k + qf + M\q ~k f 



|A-(g)| : 



(2.13) 



V±(q,k; q ,k ) 
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(k ± qf 



(k - qf + M*(q - kof 



+ (qo -* -qo) + /i 



MG 



2qk 



■In 



(fc + gf + m| 
(k - q) 2 + mj 



fl 



EL 



(2.14) 



Here /el is just a constant, and /mg is a rather complicated function of q and k. Since no singularity appears in 
/el.mg when ttid = 0, M = and q = k = /i are taken, they are considered to be sub-leading contribution and will 
be neglected in the following analysis. In the definition of V±, t he fi rst (second) term corresponds to the magnetic 
(electric) gluon contribution. We call the first integral in Eq. ( 2.13 ) as "quark-pole contribution" and the second 
integral as "antiquark-pole contribution", becaus e E-l . = q — fl (EL = q + /i) is involved in the former (latter). Near 
the Fermi surface q ~ fj,, the first integral in Eq. (2.13) gives the dominant contribution due to the small denominator 
(El + |A + | 2 )V2 = {{q _ M)2 + IA+1 2 } 1 / 2 „ 0. 

At extremely high density, the Cooper pairing is expected to take place only near the Fermi surface as we have 
discussed. In this case, we can safely neglect the antiquark-pole contribution for calculating A + . Furthermore, one 
may replace the momentum dependent vertex by the coupling constant on the Fermi surface g 2 (q, k) — > g 2 (q = 
fi,k — fi). If /i is large enough, g 2 (q — fi,k — fi) may be identified with the standard running c oup ling <? 2 (/i) = 
(87r 2 //3o)/ln(/J, 2 /Aq CD ). Thus, keeping only the magnetic interaction, one recovers the gap equation ( |2.7| ) for A + with 
g 2 — > fl 2 (/i) in the weak-coupling limit. 

At low densities, sizable diffusion of the Fermi surface occurs and the weak-cou pling approximations leading to 
Eq. (2.7) are not justified. Therefore, we need to solve the coupled gap equations Eq. (2.13) numerically. In particular, 



the replacement g (q, k) —> g < (q = fi, k = fi) is not justified when q and k are not close to /i. The contribution from 
the antiquark-pole is also not entirely negligible. 



D. Occupation number, correlation function and coherence length 

To clarify the structural change of the color superconductor from high to low densities, it is useful to examine the 
following physical quantities related to the gap function. 

(i) The quark and antiquark occupation number in each momentum state: This is a fundamental quantity character- 
izing the diffuseness of the Fermi surface. It is related to the diagonal (1-1) element of the quark propagator in the 
Nambu-Gor'kov formalism 

(^ b (t,y)^t(t,x)) supeI = km (-ilSu^ix-y)^) 



f d 3 q ■„, 
J (2^) 3 



- fa - E+iq) ) a + + - (i + ) a_ ki - PiUM 



(2.15) 



where Pg is a projection matrix to the third axis in the color space. From this expression, one can extract the quark 
and the antiquark occupation numbers as 



where the superscripts (1,2 and 3) stand for color indices. Since the third axis in the color space is chosen to break the 
color symmetry, quarks with the third color do not contribute to form Cooper pairs. When the gap is zero A± = 0, 
the system reduces to the ordinary quark matter with a sharp Fermi sphere; n]^ 2 (q) = 9([i — q) and Ti2 2 {q) = 0. 

(ii) The q-q and q-q correlation functions in the momentum space (p±{q) and in the coordinate space <p±(r): They 
reflect the internal structure of the Cooper pairs in color superconductor. These correlations are related to the 
off-diagonal (1-2) element of the quark propagator: 

(W(t,x)${t,y)) mpfa = km i[S 12 ]#(z - y) 

x°—>y°+t 

d3 Q „ia(a>-y) ( A+(g)A+(g) A_(q)A_(g) 



E±{q) 



(n 5 C)(A 2 r 2 )f?. (2.17) 



( 2 ^) 3 VV^+(<z) 2 + |A+(g)| 2 2VW+Wj 

<p±{q) is simply extracted from the above and ip±(r) is defined as the Fourier transform 

$ ± (q) = ; A±(g) p±(r) — N [ 4\ sM?) e 4<Zr , (2.18) 

2V S ±(9) 2 + |A±((?)| 2 7 (2tt) 3 

where TV is a normalization constant determined by J d 3 r|yj_|_(r)| 2 = 1. 

(iii) The coherence length £ c characterizing the typical size of a Cooper pair: It is defined simply as a root mean 
square radius of ip + (r): 

a Jrf 3 r r> + (r)| 2 / °° dfc fc 2 |rfy + (fc)/dfc| 2 
* c /d 3 r|^ + (r)| 2 JZ°dkk*\<p + (k)\ 2 ' 1 ' j 

A measure of the coherence length £ c in the weak-coupling limit is known as the Pippard length, which is given by 
£ p = (7rA + (/i)) _1 psj] . It is shown in Appendix B that the quark correlation y+(r) in the weak-coupling limit behaves 
as 

Mr-oojcxg^.e-/^). (2.20) 

In a typical type-I superconductor in metals, the Pippard length is of the semi-macroscopic order £ p ~ 10 -4 cm, 
whereas inverse Fermi momentum is of the microscopic order fcp 1 ~ 10~ 8 cm. The inverse of the Debye cutoff is in 
between the two scales lu^ 1 ~ 10~ 6 cm. Therefore there is a clear scale hierarchy, A <C wd *C kp. Because of the 
absence of the intrinsic scale wd, similar scale hierarchy in QCD at extremely high density reads fie~ c ^ 9 <C At 
lower densities, however, such scale separation becomes questionable for g is not small. 

III. NUMERICAL RESULTS 

In this section, we present numerical results of the momentum-dependent gap and the other physical quantities. 
In Sect. MA, we show solutions of the gap equations at very high de nsity . Then we discuss whether the result has 



similarity to that in the BCS superconductivity for metals [see Eq. ( 2.12 )]. What makes the color supercond ucto r 
unique is the absence of the intrinsic cutoff scale Wd- Nevertheless, we will see that similar relation as Eq. ( 2.1 2| ) 



holds under the replacement ui-o — ► A at leas t at extremely high density. We will also examine how the results of the 



gap equation in the weak-coupling limit (2/7) arc modified when the effects such as the color-electric interaction, the 



momentum dependent coupling and the antiquark-pole contribution are taken into account. 



In Sects. IIIB and IIIC, we repeat the same calculations at lower densities and show a qualitative difference from 
the weak-coupling limit. Substantial modification of the Fermi surface at low densities will be explicitly shown by 
computing the occupation number. Quark correlation in the color superconductor and the size of the Cooper pair are 
calculated in Sect. HID. The results indicate that the color superconductivity at low densities is no longer similar to 
the usual BCS-type superconductivity. 



In Sects. IIIA and IIIB, starting from the simplest gap equation in the weak-coupling limit (2.7), we will include 
the contributions of color-electric interaction, momentum-dependent coupling and antiquark pole, step by step. This 
procedure clarifies the importance of each contribution. Let us define each step below, for later convenience. 



Step 1: We solve the gap equation for A + in the weak-coupling limit Eq. ( |2.7|) where only the color-magnetic 
interaction is taken into account. The momentum independent coupling g 2 {p,p) is used in this step. 

Step 2: The Debye-screened color- electric interaction is further added to Step 1. This corresponds to solving 
Eq. ( 2.1 3| ) for A + with neglecting the antiquark-pole contribution and with making a replacement g 2 (q,k) — > 



.9 2 (A^) 

Step 3: Same as Step 2 except for the use of g 2 (q,k) instead of g 2 (n, fi). 

Step 4: The antiquark-pole contribution is also taken into account in Step 3. This gives the complete solution of our 
full coupled gap equations Eq. ( 2.1 3| ) . 



Throughout the all steps, g 2 in the Debye screening mass tod is taken to be g 2 (fj,,p) for simplicity. In Step 3 and 
Step 4, we use a phenomenological value A = 400MeV and p 2 = 1.5A 2 for g(q,k) as we have already mentioned in 
Sect. IIA. 



A. Momentum dependence of the gap at high density 



In this subsection, we solve the gap equations at chemical potential fi — 2 12 A ~ 1.6TeV. This corresponds to the 
baryon density p# = 1.1 x 10 n po with po — 0.17fm~ 3 bei ng the normal nuclear matter density. At this extremely high 
density, we expect that the full gap equation Eq. ( [2.13 ) is well approximated by its weak-coupling limit, Eq. ( [2.7] ). 
Furthermore, the analytic solution Eq. (2.9) is expected to give a fair approximation of Eq. (2.7) around the Fermi 
surface. In Fig. [j], we show the numerical results of the four steps. It is evident that all the results have a very narrow 
peak at the Fermi surface. This peak is not a singular one, but has a plateau with the width of 2A + (fc = /x), as we 
will show later. At hi gh m omentum k 3> fJ., the decrease of the gap obeys the power law. This is consistent with the 
analytic solution Eq. (2.11). Such a long tail is due to the unscreened nature of the color-magnetic interaction and is 
quite different from the behavior of the gap in the standard BCS-type superconductivity. Let us look into each step 
in more detail. 



Steps 1 and 2: 

Comparing the gap with the magnetic interaction alone (the dotted line) and the gap with both magnetic and 
electric interactions (the dashed line), we find that the electric interaction induces a large enhancement (about 10 
times) of the gap almost independently of momentum. This is already suggested by the analytic studies in Sect. IIC. 

This enhancement may be understood in a qualitative manner. In the coordinate space, the Debye-screened electric 
interaction behaves as a Yukawa potential. Such a short-range interaction can form only a loosely bound Cooper pairs 
with a large size. In fact, if one solves the gap equation with the electric interaction alone, one finds a very small gap 
compared to the one in Fig. 2(a). The situation becomes different when the magnetic and electric interactions coexist: 
Small size Cooper pairs are formed primarily by the long-range magnetic interaction. Then, even the short-range 
electric interaction becomes effective to generate further attraction between the quarks. This cooperative effect can 
qualitatively explain the reason why addition of the electric interaction enhances the gap. 

Steps 3 and 4: 

To see the effect of the momentum dependent coupling g 2 (q, fc), let us compare the dot-dashed line (where g 2 (q, k) 
is used) with the dashed line (where g 2 (n,p) is used). Since g 2 {q,k) works as a weight factor in the momentum 
integral in the gap equation, substantial difference should appear between the two cases if the contribution away from 
the Fermi surface is not negligible in the integral. However, we find no significant difference between the two cases at 
high density. This implies that the color superconductivity at high density is governed by the physics near the Fermi 
surface. 



Let us examine the last stage, Step 4, where the antiq uark- pole contribution is included. The solid line in Fig. [I] is 
A_|_ (k) obtained by solving the full coupled gap-equation ( |2.13|) for A + and A_ . Since the difference from the previous 
step is small, we can conclude that the quark-pole dominance is indeed a good approximation at very high density. 
The antiquark gap A_(fe) is shown in a separate figure (Fig. ||) together with A+(fc) (the same as the solid line in 
Fig. 1). Unlike A+(fc), the antiquark gap A_(fc) is a smoothly decreasing function of the momentum. Although A_ 
is not small compared to A+, it does not imply that the sizable antiquark Cooper pairs exist in the system because 
the number of antiquarks are much smaller than quarks as we will show in Sect. HID. 
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FIG. 1. The momentum dependence of the gap A+ from Step 1 to Step 4 at high density, ^ = 2 12 A. The dotted line is 
a result of the magnetic interaction alone (Stepl). The dashed line is the gap with both magnetic and electric-interactions 
included (Step2). The dot-dashed line is the gap in which the momentum dependent coupling g 2 (q, k) is used (Step 3) in stead 
of the momentum independent coupling g 2 (fi, fi) (Step 2). The solid line is the gap as a solution of the full gap equation (2.13) 
with g 2 (q,k) and the antiquark pole (Step 4). 

Finally, Fig. |^ shows a comparison between the momentum dependence of the numerical solution A + in St ep 4 (the 
solid line which is the same as the solid line in Fig. ||) with the analytic solution near the Fermi surface (2.9) (the 
dashed line). The peak height of the analytic solution is adjusted so that the solid and dashed lines coincide at k = fj,. 
In the vicinity of the Fermi surface (Fig. 3(a)), the momentum dependence of the analytic solution agrees with the 
numerical solution reasonably well. Note that the peak is not a singular one as we mentioned before. On the other 
hand, as is shown in Fig. 3(b) the deviation becomes considerable away from the Fermi surface. In particular, the 
asymmetry with respect to the Fermi surface k — /i represents the deviation from weak coupling analytic solution. 
Note that the coincidence at k ~ is accidental because the analytic result ( |2.9| ) is valid only near the Fermi surface. 

In this subsection, we have solved the momentum-dependent gap equations at high density /i = 2 12 A. The char- 
acteristic features of the momentum-dependent gap are (i) there is a sharp peak at the Fermi surface, and (ii) the 
gap decays rapidly but is nonzero for mo mentu m far away from the Fermi surface. The property (i) is similar to the 
standard BCS superconductivity [see Eq. (2.12)] but (ii) is not, due to the absence of intrinsic ultraviolet Debye-cutoff 
of the gluonic interaction in QCD. As for the magnitude of the gap at high density, the color-electric interaction 
enhances the gap considerably. The effects of the momentum-dependent coupling g 2 (q,k) and the antiquark pairing 
are shown to be not important at high density. 
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FIG. 2. A+ and A_ as the solutions to the full coupled gap-equation in Step 4 at high density, fi — 2 12 A. 
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FIG. 3. Comparison of the numerical solution of A+ in Step 4 (the solid line) and the analytic solution ( [2,S| ) (the dashed 
line) at high density, /i = 2 12 A. The latter is normalized to the former at k — fi. (a) near the Fermi surface, (b) global behavior. 
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B. Momentum dependence of the gap at low density 



In this subsection, we carry out the same calculations as in Sect. IIIA at lower density \i = 2A = 800MeV (which 
corresponds to the baryon density pe = 13.2po)- In Fig- @, we show the numerical results of the four steps. One of 
the main differences from the high density case is the absence of a sharp peak. The width of the peak in Fig. ^ in 
terms of k/fi is much wider than that in Fig. [I]. This large modification of the gap function will be attributed to the 
larger coupling at low density, because the Cooper pairs far from the Fermi surface are easily formed. Therefore, this 
implies that the physics of superconductivity is no longer limited on the Fermi surface. We confirm this below by 
more detailed analysis of each steps. 
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FIG. 4. Same as Fig. |l| except that the density is lower, fi = 2A = 800MeV. 

Steps 1 and 2: 

Comparing the dotted line (with the magnetic interaction alone) and the dashed line (magnetic and electric inter- 
actions) in Fig. ^, we again find that the addition of the color-electric interaction induces a large enhancement of the 
gap. 

Steps 3 and 4: 

The result of Step 3 (where g 2 (q,k) is used) is shown as the dot-dashed line, which should be compared with the 
previous step (the dashed line) where g 2 (/i,/x) is used. In contrast to the high density case in Fig. |], the momentum 
dependent coupling makes the magnitude of the gap considerably smaller. This is one of the evidences that the 
momentum integration is not any more dominated by the small region near the Fermi surface. The reason why the 
suppression of the gap takes place instead of the enhancement is understood as follows. Consider the gap at the Fermi 
surface A + (fc = /j). In the momentum integration of the gap equation Eq. ( 2.13] ) at fc = jj,, g 2 (q,k — (i) < g 2 ((A, fi) for 
jj, < q, while g 2 (q,k = /x) = g 2 (/i,//) for /i > q. Therefore, the momentum-dependent coupling always acts to reduce 
the magnitude of the integral as compared to the momentum-independent one. 

Now , let us include the antiquark excitation in the gap equation (Step 4). The solution of our full gap equation 
(243) is shown by the solid line. The antiquark-pole contribution enhances the gap by 10% . The difference between 
the solid and dot-dashed lines is not huge but the absolute magnitude of the enhancement is larger than that in the 
extremely high density case shown in Fig. 1. The reason for this enhancement is clear: As one decreases density, the 
chemical potential becomes small, which plays a role of the intrinsic energy gap for antiquark excitation. Therefore, 
the antiquark-pole contribution becomes non-negligible. This also suggests that the more low-momentum antiquarks 
are present at low density which will be confirmed in Sect. HID by computing the occupation number of antiquarks. 

Finally, in Fig. ||, we compare the quark-gap A + (the solid line) and the antiquark gap A_ (the dashed line). As 
in the case of high density, the antiquark gap A_ is a smooth function of k and is not small compared to A + . 



In this subsection, we found that the momentum dependence of the gap at low density is quite different from th at a t 
very high density. The sharp peak at the Fermi surface disappears. The gap equation in the weak-coupling limit (2.7) 



is no longer a good approximation and all the contributions neglected in the weak-coupling limit are not negligible. 
All these results allow us to conclude that the color superconductivity at low density is not a phenomenon just around 
the Fermi surface. In subsections HID and HIE below, we will strengthen this picture by looking at other quantities 
such as the quark occupation numbers and the size of the Cooper pair. 
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FIG. 5. Same as Fig. || except that the density is lower, fi = 2A = 800MeV. 



C. Density dependence of the gap 

Let us now examine how the gap at high density with a sharp peak changes into the gap at low density with only 
a broad bump. In Fig. ||(a), we show the gap A + (fc) as a solution of the full gap equation (Step 4) for a wide range 
of densities. Since the actual position of the Fermi surface moves as we vary the density, we use k/fi as a horizontal 
axis in the figure in order to show the change of the global behavior. The figure shows that the sharp peak at high 
density gradually gets broadened and simultaneously the magnitude of the gap increases as we decrease the density. 

In Fig. |^(b), the gap at the Fermi surface A + (/i) is shown as a function of the chemical potential. It decreases 
monotonically as \i increases, but turns into an increase for fj, > 10 6 MeV. The analytic solution is also shown in 
Fig. |^(b). The magnitude of the analytic solution is normalized to the numerical solution at the highest density 
/i = 2 12 A ~ 1.6 x 10 6 MeV. At high density, /i-dependence of the numerical result is in good agreement with the 
analytic form which has a parametric dependence A + (/i) oc g~ 5 /j,exp(— 37r 2 /v2g) with g 2 = g 2 (/i,^). On the other 
hand, the difference of the two curves at low density implies the failure of the weak-coupling approximation. As we 
have seen before, the use of g 2 (q,k) and the antiquark-pole has non- negligible effects on the gap in the low density 
regime. 
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FIG. 6. (a) A+(fc) as a function of k/fj. for various densities = 2 n A with n = 1, 2, 3, 12. All the calculations are done with 
the momentum dependent vertex g 2 (q, k) and with the antiquark-pole contribution, (b) Chemical potential dependence of the 
gap A+(fc = fi) in the full calculation compared with the analytic result which is normalized at the highest density /i = 2 12 A. 



D. Occupation numbers 



So far we have seen that the weak-coupling picture of the color superconductivity is modified at low densities. 
In order to see how the Fermi surface is diffused by the Cooper pairing, let us evaluate the occupation numbers of 
quarks and antiquarks. Since w e hav e solved the full gap equation (2.13) for A±(fc), we can immediately obtain the 
occupation numbers using Eq. (2.16). 

In Fig. |?], the occupation numbers of quarks [Fig. 0(a)] and antiquarks [Fig. |7|(b)] are shown at high density p, = 2 12 A 
and at lower densities /i = 2A, 2 2 A. One finds that the quark occupation numbers is almost a step function at high 
density, while it is smeared out for a wide region of momentum at low densities. Th e diff useness of the Fermi surface 
is found to be of the order of A+(/z), which is consistent with the definition in Eq. (2.16). 

Figure @(b) implies that a small amount of the antiquarks also participate in the color superconductivity. Although 
the antiquark gap A_ is of the same order of the quark gap A+, the antiquark occupation number is generally 
suppressed due to the large energy denominator {(k + /i) 2 + IA-1 2 } 1 / 2 . As one decreases the density, however, such 
suppression is relatively weakened and the magnitude of the antiquark occupation number increases, as can be seen 
from Figure 0(b) . 




FIG. 7. Occupation numbers of quarks nV 2 (a) and antiquarks n^ 2 (b) for several different densities. 



E. Correlation function and coherence length 



One of the advantages of treating the momentum dependent gap is that we are able to calculate the correlation 
function which physically corresponds to the "wavefunction" of the Cooper pair. Such correlations have been first 
studied in Re f. [ 15f| in the context of the color superconductivity. 

Using Eq. (|2.18f ), we calculate the correlation functions of quarks and antiquarks, <f±(k). The results at various 
densities are shown in Fig. pi For quarks [Fig. ||(a)], the correlation function at very high density has a sharp peak 
at the Fermi surface but it becomes broader as we decrease the density. This is of course due to the broadening of 
the gap which we found in Fig. |6|. 

For antiquarks [Fig. ||(b)] , the correlation is much weaker than that of quarks and is a smoothly decreasing function 
of k. Also, the magnitude of the correlation increases as we decrease the density. Since A_(/c) <C /i holds for 
densities considered in this paper, the above features can be simply understood by an approximate relation, ip~(k) ~ 
A_(*)/2(fc + A*). 
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FIG. 8. The correlation functions in the momentum space at several different densities for quarks (a) and for antiquarks (b) . 



We can compute the size of a Cooper pair from the correlation function. Figure |^(a) shows the co heren ce length £ c 
of a quark Cooper pair defined as the root mean square radius of the correlation function [see Eq. ( 2.1S| )1. The size 
of a Cooper pair becomes smaller as we go to lower densities. This tendency is understood by the behavior of the 
Pippard length £ p = 1 hx A + (/i) (which gives a rough estimate of the coherence length) together with the behavior of 
A + (/x) shown in Fig. ||(b). 

Also, the Cooper pair becomes smaller as density increases beyond /i = 2 12 A. However, it does not necessarily 
imply the existence of tightly bound Cooper pairs. In fact, the size of a Cooper pair makes sense only in comparison 
to the typical length scale of the system, namely the averaged inter-quark distance d q defined as ^ 

As we go to higher densities, the ratio £ c /d q increases monotonically as shown in Fig. ^|(b). Namely loosely bound 
Cooper pairs similar to the BCS superconductivity in metals are formed at extremely high densities. 

At the lowest density in Fig. g, the size of the Cooper pair is less than 4 fm and the ratio £, c /d q is less than 
10. The transition from £ c /d q ^ 1 to £ c /d q ~ 1 as fj, decreases is analogous to the transition from the BCS-type 



2 This is a result of free quarks. To obtain accurate density p and d q , we have to include contributions from interaction. 
However, such correction is suppressed by the power of A+ //i. 



superconductor to the so-called "strong coupling" superconductor. The weak-coupling BCS superconductivity may 
smoothly change into the Bose-Einstein condensation (BEC) of tightly bound Cooper pairs as the coupling strength 
increases p6| . Our result here suggests that the quark matter possibly realized in the core of neutron stars may be 
rather like the BEC of tightly bound Cooper pairs. 

For better understanding of the internal structure of the quark Cooper pair, let us consider the correlation function in 
the coordinate space. Fig. shows the spatial correlation of a Cooper pair at various chemical potentials normalized 
as J d 3 r\if + (r)\ 2 = 1. As is expected, the density dependence of the quark correlation in the coordinate space is 
opposite to that in the momentum space. At high density, most of the quarks participating in forming a Cooper pair 
have the Fermi momentum kp = fi giving a sharp peak in the momentum space correlation. In the coordinate space, 
this corresponds to an oscillatory distribution with a wavelength A = without much structur e nea r the origin. 
(The oscillation is also evident from the factor sin(/Ltr) in the approximate correlation function Eq. (2.20) discussed in 
Sect. HE.) At lower densities, accumulation of the correlation near the origin in the coordinate space is much more 
prominent in Fig. [l0| This implies a localized Cooper pair composed of quarks with various momentum. This is also 
seen by the broad momentum correlation in Fig. pKa). 





IV. SUMMARY AND DISCUSSION 



In this paper, we have studied the spatial-momentum dependence of a superconducting gap and the structure of the 
Cooper pairs in two-flavor color superconductivity, using a single model for a very wide region of density. Nontrivial 
momentum dependence of the gap manifests itself at low densities, where relatively large QCD coupling allows the 
Cooper pairing to take place in a wide region around the Fermi surface. Our results imply that the quark matter 
which might exist in the core of neutron stars or in the quark stars could be rather different from that expected from 
the weak-coupling BCS picture. 

Following is the summary of what we have discussed in this paper. 

(1) At high density, the weak-coupling gap equation with the electric and magnetic gluons is a good approximation. 
The momentum dependent quark-gluon vertex and the coupling to the antiquark pairing do not change the weak- 
coupling result. The gap A_|_(fc) has a peak around the Fermi surface with a width ~ 2A+(/x) and decreases 
rapidly as k goes away from the Fermi surface. This is consistent with the analytic solution of the gap equation 
in the weak-coupling limit. All the results indicate that the Cooper pairing takes place only in a small region 
around the Fermi surface. 

(2) At lower density, the weak-coupling gap equation is no longer a good approximation. The momentum-dependent 
vertex and the coupling to the antiquark pairing have non-negligible effects. The sharp peak at the Fermi surface 
disappears. These imply that a large number of Cooper pairs away from the Fermi surface participate in the 
color superconductivity. This was confirmed by the quark occupation numbers and the quark-quark correlation. 
Therefore, color superconductivity at low density is not a phenomenon in the vicinity of the Fermi surface, but 
is a phenomenon with large modification of the Fermi surface. 

(3) The qualitative change of color superconductivity from high density to low density can be explicitly seen by the 
ratio of the size of the quark Cooper pair to the averaged inter-quark distance. At high density, the ratio is 
very large (about 10 5 at /j, <~ 10 6 MeV) which is consistent with the standard BCS-type superconductivity. At 
lower densities, however, the ratio becomes small (about 10 at /i = 800MeV). This situation is rather similar to 
the "strong coupling" superconductor which could be described by a Bose-Einstein condensate of tightly bound 
Cooper pairs. 

There are several future problems. 

Firstly, there are still several corrections to our "full" gap equation. They include the full hard-dense-loop corrections 
and the Meissner effect in the gluon propagator, and also the use of the plasmino dispersion in the diagonal self-energy 
of the quark propagator. The latter effect induces only a sub- leading change on the gap at the Fermi surface at high 
density p^ps] , but may have non-negligible effects at lower densities. For treating both the diagonal and off-diagonal 
parts of the quark self-energy, the standard Eliashberg formalism must be used |^,^3|. If we include all the above 
effects, the gap equation with frequency and momentum as independent variables should be solved J3(J. 

Secondly, we assumed that there is no significant vacuum effects in the present paper. If the density is close to the 
critical density of chiral symmetry brea king , one must consider the interplay between the quark-antiquark condensate 
and the quark-quark condensate |pH , p^^lf . This also requires us to treat the diagonal and off-diagonal components 
of the fermion self-energy simultaneously, which corresponds to the Hartree-Fock-Bogoliubov theory in many-body 
problem [||. 

Thirdly, we have taken the Landau gauge £ = in this paper. As we discussed in Sect. IID and in Appendix A, the 
^-dependence in the gap equation remains as far as we properly treat the momentum dependence in the gap equation. 
Although it is a sub-leading contribution relative to the magnetic and magnetic interactions, it is desirable to develop 
a gauge invariant formalism especially when one treats the low density region. 

Lastly, it will be very intriguing to find a new way of describing color superconductivity at low densities. The small 
size Cooper pairs with coherence length comparable to the inter-quark distance suggests that BEC description may 
be useful as in the analogous example in condensed matter physics. An analysis along this line has been discussed in 
Ref. H using the linear sigma model for the diquark field. 
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APPENDIX A: GAUGE DEPENDENCE OF A 



In this Appendix, we discuss the gauge parameter dependence of the gap equation for A + at high density. We write 
the gauge dependent contribution to the gap equation for A+(&) as if(fc). Neglecting the antiquark pole, it becomes 



24tt 2 



A 1 



1 ■>{-<>) - A+(g) 



K^{q, k; cost 



(Al) 



where q ■ k — cos 9 and K$ is defined as 



K^{q,k;x) = £ x — 
q 2 



2qk 



(q — k cos 9)(k — qcosff) 1 — cos 9 



k 2 — 2qk cos 9 \ q 2 + k 2 — 2qk cos 9 



(A2) 



I^(k) should be added to the right-hand side of Eq. ( |2.13| ) when £ ^ 0. 

Now, if we set q = k = fi before the angular integration, the kernel vanishes for all cosf?, namely K^(p, /i; cos#) = 0. 
This is the standard argument that the gauge parameter dependence does not appear in the gap equation at extremely 
high density. However, if we first integrate over cos 9, and then take the limit q, k — > /i, the result is nonzero. This 
can be seen explicitly by carrying out the angular integral of and by writing the result in terms of a variable 
Y = (k/q + q/k)/2 > 1 : 



J d(cos9) K%(q, k; cosO) = £ 



Y- 1 



In 



Y + l 



Y - 1 



(A3) 



If we take the limit Y — * 1 in the right-hand side, the first term disappears, but the second term survives and gives a 
gauge-dependent contribution — £. One can trace back the origin of this situation by introducing a small regulator e 
to the collinear region of the integral cos 9 ~ 1. Then, the integral above becomes 



F{e,Y) 



e?(cos 9)K^(q, k; cos ( 



Y 



In 



Y - 1 



Y - 1 Y + l 
Y-l+e 2 



(A4) 



The first term disappears irrespective of the order of Y 
when one takes Y — ► 1 before e — > 0. Namely, 



1 and e — ► +0. However, the second term vanishes only 



lim lim Fit, Y) = 0, 

lim lim Fit, Y) = £ . 
y^i e^o 



(A5) 
(A6) 



This non-commutability of the two limits arises from the fact that F(t, Y) does not approach to F(0, Y) uniformly 
in the region Y > 1. 

Therefore, if we integrate over the angular variable cos 9 exactly, the gap equation for A + (fc) at high density has 
an extra gauge-dependent contribution from the quark-pole; 
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dq 



^/E + (q) 2 + A + (q) 2 k 



4qk 



{k + qf 



(k - qf 



(A7) 



Compared with the leading magnetic contribution with the kernel of logarithmic singularity, the above contribution 
is considered to be a sub- leading effect. 



APPENDIX B: PIPPARD LENGTH IN THE WEAK-COUPLING LIMIT 



In this Appendix, we derive the Pippard length from the corre lation function in the weak-coupling region. We use 
the two-point correlation function of quarks defined in Eq. (2.1 8|) 



<p+(q) 



A+(g) 



Since <p+{q) is a function of q(> 0) only, its Fourier transformation is given by 



(Bl) 



N < • _2. 



<P+{r) = 2~2 J d 11 3o(qr)<p+(q) 



2i 2 r Jo ^(,-rt2 + A+(«) ! 

where jo(a;) = x _1 sinx is the 0-th order spherical Bessel function. 

In the weak-coupling limit, the integral is dominated at q ~ fj,. Then, by replacing A+(g) by A+(/x) and performing 
the g integration approximately, one finds 

r \ « f , . M MA + (M) 



2 ^ V(9-a*) 2 + a+( m ) 

A^/x 3 sin(/xr) A + (/i) 



7r 2 /ur /i 



X (A + (/i)r), (B3) 



where Kq is the 0-th modified Bessel function of the second kind. Using its asymptotic form 

,oo nr r(n + i/2) 1 



olZj ~ \Yz e 2^r(-n + l/2)n!(2z)™' 

f n=0 v / y \ / 



we obtain 



2tt 3 (ttr) 3 / 2 

where £ p = (7rA + (^i)) _1 is the Pippard length. 
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